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Background 


Issues; 

• COPVs can be at risk for catastrophic failure 

- Risk of insidious burst-before-leak (BBL) stress rupture^ (SR) failure of 
carbon-epoxy (C/Ep) COPVs during mid to late life 

- Risk of lower burst strength of C/Ep COPVs subjected to impact damage 

• Issues with manufacturing defects and inspectabilty of COPVs on 
NASA spacecraft (ISS, deep space) 

• Lack of quantitative NDE is causing problems in current and future 
spacecraft applications 

- Must increase safety factor or accept more risk 

- Thinner liners are driving need for better flaw detection in liner and 
overwrap 

^ SR defined by AIAA Aerospace Pressure Vessels Standards Working Group as “the 
minimum time during which the composite maintains structural integrity considering 
the combined effects of stress level(s), time at stress level(s), and associated 
environment. ” 
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STRAIN <%) 


Strain vs. Tl 


Behavior During Creep 



Classical Case 



t, Time 


distinct tertiary creep phase 

(ductility observed before rupture) 



BBL rupture of 


aCOPV 


The probiem with advanced fibers such as carbon or Keviar® is that no 
ductiiity is observed before rupture during tertiary creep, so the 
stress rupture occurs with iittie or no advance warning 
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Effect of Fiber Choice on Stress Rupture 



C/Ep COPVs are susceptible to stress rupture, although to a lesser extent 

than glass or Kevlar® fiber composites 
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Characteristic lifetimes of graphite, Kevlar® and glass-reinforced composites at 
different percentages of the ultimate strength. Each symbol represents the 
median life (50%) under sustained loads as percentage of the ultimate strength 

of the material § 


§ 


J. T. Shaffer, ‘‘Stress Rupture of Carbon Fiber Composite Materials'', in Proc. 18^^ Inti. SAMPE Technical Conference, p. 613 (1986). 
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COPVs on ISS 


• Presently have 1 7 high pressure COPVs on ISS (most are C/Ep) 

- Up to seven additional COPVs are planned and under development 

• Long term reliability risk levels are 1 0 ® or lower except for NTA and 
SpaceDRUMS COPVs, which have risk levels of 10'^ to 10'^ § 

- Reliability much lower if C/Ep overwrap sustains impact damage 


Subsystem 

No. 

Shape 

Size. in. 

Commodity 

Materials 

Supplier 

FOS 

MEOP psi 

Liner 

Wrap 

ECLSS.ACS 

HPGT 

4 

Sphere 

37.89 

Ox>gen. 

Nitrogen 

301 SS 

IM.7W 

GD 

2.0 

5000 

ECLSS/MCA 

1 

Cylmder 

7.22 Lx 3.55 D 

Calibrated air 

A1 

S-Glass 

SCI 

3.4 

3000 

TCSOTA 

2 

Cylinder 

45 Lx 19.7 D 

Nitrogen 

A1 

T-1000 

GD 

2.52 

3000 

E\^SAFER 

3 

Cylmder 

9Lx6D 

Nitrogen 

SS 

T-1000 

ARDE 

3.0 

10.000 

Emiromnents P 
CU 

2 

Sphere 

15.37 

Xenon 

301 SS 

T-1000 

ARDE 

4.17 

3000 

Payloads. 

SpaceDRUMS 

5 

Cylmder 

17.1 L X 8.5 D 

Argon 

A1 

T-1000 

GD 

2.28- 

2350 

Pa>toads. 

VCAM* 

1 

Cylmder 

8.1 Lx 3.68 D 

Hebiim 

A1 

Gr ep-2150 

Carleton 

3.4 

1985 

AMS-02* 

2 

Sphere 

12.4; 15.8 

CarbiMi dioxide. 
Xenon 

301 SS 

T-1000 

ARDE 

3.05-4.4 

1440-2900 

ECLSS&TCS/N 

ORS^* 

0 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

TBD 

CPAA^*^ 

4 

Cylmder 

4.04 Dx 9.6 L 

Carbon dioxide 

A1 

Gr E-Glass 

Carleton 

4.67 

4500 


*The VC AM and AMS systems have not been manifested. 



**The NORS system is still imder de\’elopment. 

^*^The CIPAA system is transported to and from the ISS with each Shuttle mission. The very high FOS indicates a very 
low risk of rupture. 


E. Y. Robinson, R. Kohli, ‘'Preliminary Stress Rupture Risk Assessment for Graphite/Epoxy Composite Overwrapped Pressure Vessels on 
the International Space Station ”, Aerospace Report No. ATR-2009(5298)-6, Sept. 30, 2009. 
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Goals 


• Develop quantitative AE procedures specific to C/Ep 
overwraps, but which also have utility for monitoring 
damage accumulation in composites in general 

• Lay groundwork for establishing critical thresholds for 
accumulated damage in composite components, such 
as COPVs, so that precautionary or preemptive 
engineering steps can be implemented to minimize or 
obviate the risk of catastrophic failure 

- Felicity ratio (FR), coupled with fast Fourier transform (FFT) 
frequency analysis shows promise as an analytical pass/fail 
criterion 

• Would fail COPVs at a critical FR (FR*) below 1 .0, indicative of 
severe accumulated damage 

• Could also fail COPVs at a known levels cumulative of fiber 
breakage or matrix cracking 



Load control and AE data acquisition system (DAGS) consisted of: 

• Instron® 5569 Series Electromechanical Test Instrument (left) 

• DigitalWave Corp. FM-1 8-channel DAGS (lower right) 

• AE and tensile test GPU controllers (upper right) 
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Experimental 


AE Sensors: Each channel (4 used) was connected to a DWG PA-0, 0 dB Gain 
preamplifier, and then to a broadband high fidelity B1080 piezoelectric sensor 
with a frequency range 1 kHz to 1 .5 MHz. Sensors were mounted on 
cardboard-tabbed C/Ep tow specimens (8-in. gage length) using Lord Corp. 
AE-10 acrylic adhesive 



Experimental 


Tabbing: shear strength of epoxy and bonded grip length important variables^ 



= rw2P“ 


where: 

^min = minimum required bonded tab length, mm fin.]; 

= ultimate tensile strength of coupon material. MPa 
[psij: 

h = coupon thickness, mm [in.]; and 

= ultimate shear strength of adhesive, coupon mate- 
rial. or tab material (whichever is lowest), MPa 
fpsi]. 




s 

ASTM D 2343, Test Method for Tensile Properties of Glass Fiber Strands, Yarns, and Rovings Used in Reinforced Plastics, American 
Society for Testing and Materials, West Conshohocken, PA (2008). 

ASTM D 3039, Test Method for Determining Tensile Properties of Polymer Matrix Composite Materials, American Society for Testing 
and Materials, West Conshohocken, PA (2007). 


Felicity Ratio Analysis 

(IM7 & T1 000 composite tow) 


Experimental 


• For purposes of quick turn-around, an intermittent ioad hoid (ILH) 
stress schedule was used (red data) 



• iLH profiie is based on the pressure tank examination procedure 
described in ASTM E 1 067 § 

S 

ASTM E1067, Practice for Acoustic Emission Examination of Eiberglass Reinforced Plastic Resin (ERP) Tanks A^ess els, American Society for 
Testing and Materials, West Conshohocken, PA, 19428- (2011). 


FELICITY RATIO 


Results & Discussion 

composite tow 


Linear decrease in FR with load noted for T1 000 (R^ > 0.9) and IM-7 (R^ > 0.99) 

C/Ep, similar to the behavior noted for Kevlar 49-epoxy K/Ep 

. : 
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• For same material and averaging method, the slope of least squares fit is 

indicative of damage tolerance 

- Flatter slopes correspond to good damage tolerance (in-character behavior) 

- Steep slopes correspond to low damage tolerance (out-of-character behavior) 

• Kaiser effect violated at FR<1 ^ onset of severe accumulated damage 

• C/Ep produced more AE than K/Ep (but AE less energetic on average) 12 


Results & Discussion 

composite tow 



Formation of characteristic damage state very evident at ioad ratio 
(LR) < 0.6 for IM-7 


« MesjnCnmfir 


- L*«f*roAlr(KM<l| — fot4»Cw^l 



- In quasi-isotropic composite lay-ups, for example, characteristic damage state 
formation thought to involve predominant matrix cracking 

- For uniaxial tow, FFTs revealed the characteristic damage state formation involves 
mixed mode failure mechanism (cooperative matrix cracking, fiber/matrix 
debonding, fiber pull-out, fiber breakage) 


Results & Discussion 

composite tow 


Date 

Material 
& Spool # 

Filteri 

F@ 

FR=1 (lb,) 

F 

max 

(lb,) 

FR=\ 

(ksi) 

®max 

(ksi) 

FR^ 

Failure 

Type2 


83109 

IM7#95 

32% 

135 

210 

342 

532 

0.95 

XGB 


90109 

1M7#95 

27% 

151 

234 

383 

591 

0.945 

XGM 


90809 

IM7#95 

58% 

171 

210 

433 

530 

0.971 

XGM 


111009 

IM7#117 

9%3 

193 

252 

488 

637 

0.961 

XGM 


32610 

IM7#61 

19% 

183 

228 

464 

578 

0.97 

XGM 


82509 

TIOOO #74 

32% 

240 

355 

658 

972 

0.972 

XGT 


82609 

TIOOO #74 

46% 

231 

369 

633 

1010 

0.953 

XGT 


82809 

TIOOO #74 

37% 

226 

362 

UTS 

992 

0.977 

XGT 


112409 

TIOOO #155 

4%3 

181 

379 



0.945 

SGM 



1037 


112509 

TIOOO #74 

6%3 

206 

325 

WaO / aW /O 

scatter 

890 

0.966 

LGM 


40910 

TIOOO #155 

6%3 

181 

374 

493 

1024 

0.95 

XGM 


Mean 

IM7 

29% 

167 

111 

422 

/' 57S''\ 

/ 0.959'\ 


FR* 


Std. Dev. 

18% 

24 

18 

60 

/ \ 
; 45 ; 

\ 0.012 ; 

1.2-1 .4% 

Mean 

TIOOO 

22% 

211 

361 

577 

\ 988 /' 

\ 0.961 /' 

scatter 


Std. Dev. 

18% 

26 

19 

71 

T~ 

53./ 

7 — 

0.013. 




• Let FR^ = extrapolated FR at rupture predicted by the least squares fit 

• FR" behaves like a universal parameter that varies less than the UTS 

^ Data filter reflects percentage of events removed from the raw AE data 

2 Failure abbreviations per ASTM D 3039, Test Method for Determining Tensiie Properties of Poiymer Matrix Composite Materiais, 

American Society for Testing and Materials, West Conshohocken, PA (2007) 

3 Improved tabbing method 



Results & Discussion 

composite tow 

Source location of FR events show they occur at or 

near locus of failure 


• IM7 03261 0 specimen had intact tow 
between and 0 (lower tab) and 0.11 5 m 
0.17 and 0.20 m (upper tab) 

• Tow region between 0.1 15 and 0.17 m 
obliterated (explosive failure) 

• Most FR events were source located in 
the missing region that failed 
explosively 


IM7_032610 FR Event 
Locations 
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Conclusions 

composite tow 

• Consistent FFT values noted for T 1 000 and IM-7 

• Suggests that the FR can be used as an analytical 
PASS/FAIL criterion for C/Ep composite materials 

• Precedent: ASTM suggests using FR < 0.95 as failure 
criteria in fiberglass reinforced pressure vessels § 

- Experimental C/Ep failure criteria from strand tests 

» IM7: FR< 0.959 

» T1000: FR< 0.961 

- Also can use counts and number of hits above high energy threshold 

• Opens up possibility that C/Ep composite materials can 
be subjected to ILH profiles to assess in- or out-of-family 
response 

- Need to verify that test specimens or articles with low initial FR, or steep 
‘FR vs. load’ slopes in fact fail prematurely, or in the case in COPVs, fail 
at lower burst pressure 

^ ASTM E1067, Practice for Acoustic Emission Examination of Eiberglass Reinforced Plastic Resin (ERP) TanksP/essels, American Society for 
Testing and Materials, West Conshohocken, PA, 19428 (2011). 



Waveform and FFT Analysis 

(IM7 & T 1 000 composite tow) 
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Results & Discussion 

Waveform and FFT Analysis 

AE frequency ranges have been correlated with 
micromechanical damage mechanisms in C/Ep§ 




Peak Frequency 


■ Fiber Breakage- 

■ Debonding 

■ Fiber Pullout — | 

■ Matrix Cracking 



De Groot, P., P. Wijnen, and R. Janssen, “Real-time Frequency Determination of Acoustic Emission for Different Fracture Mechanisms in 
Carbon/Epoxy Composites,” Composites Sci. TechnoL, 55, pp. 405-421 (1995). 

Dzenzis, Y. A., and J. Qian, “Analysis of Microdamage Evolution Histories in Composites,” Int. J. Solids and Structures, 38 , 1831-1854 (2001). 


Amplitude, Volts 


Results & Discussion 

Waveform and FFT Analysis 

FFT (unfiltered) showing concerted failure using De Groot’s 

frequency ranges 



FFT FREQUENCY DISTRIBUTION 

TIOOO Spool 74 tested 9/9/09, Y=14.8 cm (2/5 from S3 to S4) 
N=2597, E=3.39 FAC-4 
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Results ^ discussion 

Waveform and FFT Analysis 


• In general, three different waveforms were observed for C/Ep 

1. Matrix Cracking 2. Fiber Breakage 
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Results Si Discussion 

Waveform and FFT Analysis 



Three different waveforms were observed for C/Ep 

3. Concerted, mixed mode failure 
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Results & Discussion 

Waveform and FFT Analysis 


• IM-7 early vs. late life events 


Early 

Life 



Late 

Life 



LR < 0.5 

first ILH ramp up 


LR > 0.8 

last ILH ramp up 


Notice change from ordered (early) to unordered peaks (late life) 


22 


Results & Discussion 

Waveform and FFT Analysis 


• High frequency peaks shifted downward with increasing load ratio 

731 kHz ^ 728 kHz ^ 685 kHz ^ 640 kHz 

• Attributed to increasing accumulated damage, hence lower 
modulus, causing siower stress wave propagation 



Results & Discussion 

Waveform and FFT Analysis 



• The FFTs of IM-7 and T 1 000 Felicity ratio events (first ten events) 
were then compared to see if they had a characteristic damage mode, 
or if the damage mode changed with ioad 



• Fiber breakage dominates FR events 

- Otherwise FR events involve concerted failure for both types of C/Ep 

• Some differences, but same overail trend noted for T1 000 & IM-7: 

300-1000 kHz > 90-190 kHz > 190-300 kHz 
(fiber breakage > matrix cracking > debonding/puii-out) 
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Data Reduction 
Enhancements 
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Results & Discussion 

Data Reduction Enhancements 


• Felicity Ratio Analysis Tool (FRAT) written 

- Automate AE data reduction 

- Optimizes best fit using least squares or bisquare fitting 








Results & Discussion 

Data Reduction Enhancements 



• Filtering Criteria 

- -1 arrival channel events were plotted against the load profile 

- -1 events primarily occurred on loading & upper load holds 

- No grounds for rejecting lower energy events 


Results & Discussion 

Data Reduction Enhancements 


Linear Least Squares 
(LLS) 

- Gives outliers too much leverage 
over fit. 

- Must manually remove outliers 

- Minimizes square of residuals 

Bisquare weighting 

- Very similar process to LLS 

- Weights residuals of each point 
and down weights points of high 
leverage. 

• Automatically remove outliers 

- Minimizes weighted square of 
residuals 
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Results & Discussion 

Data Reduction Enhancements 



Can optimize the linear fit for different methods for determining 
the onset of significant AE used to calculate the FR: 



best linear fits 


■ , ! .j. - 

: 1 



■ "t 't' “■ ■ 

■ t • ; 



r* 

r ! 

■ t ■: 

: 


P' 

CD* 1 




n 




FR^ 




n% 


n-bar % 


29 


Critical Felicity Ratio. FR* 


Results & Discussion 

Data Reduction Enhancements 


• Comparison of IM7 and T 1 000 tow showing the variation of the 
(coefficient of determination) and FR* with n: 



0 10 20 30 40 so 



10 

09 

05 
07 

06 
05 
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n-bar 


0 10 20 30 40 so 



Note: iower values of R^ for T 1 000 as compared to IM7 
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Results & Discussion 

Data Reduction Enhancements 


• Analytical identification of the ‘knee’ § in the AE events vs. time curve 
using an exponentially weighted moving average (EWMA) method: 





Results & Discussion 

Data Reduction Enhancements 


• The ‘knee’ is not analytically defined in ASTM E2478§ necessitating 
development of an exponentially weighted moving average (EWMA) 
method: 
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Zj = IXi + (1 — ^)Zj_i 
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S ^ 

Limits = T±L— I- r[l-Cl-l)^'] 

2 — A 


standard deviation: 

" 0-iK(2) 


ASTM E2478, Standard Practice for Determining Damage-Based Design Stress for Glass Fiber Reinforced Plastic ( GFRP) Materials Using Acoustic Emission, 
American Society for Testing and Materials, West Conshohocken, PA, 19428 (201 1). 
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Results & Discussion 

Data Reduction Enhancements 


• EWMA method found to yield better (more linear) results than other 
methods for T 1 000 tow: 



better 
linear fit 
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Results & 

Data Reduction Enhancements 


Effect of data richness on FR trend analysis linearity showing that 
more sensitive AE DACS setting produce better linearity: 


1.0 


I 0.9 

16 

c 


I 0.8 
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For example, 30 dB sensitivity is better than 50-60 dB sensitivity 
for FR analysis or quantification of fiber breakage events 


Application to 

Composite Overwrapped Pressure Vessels 

(COPVs) 


Results & Discussion 

COPVs 


A 6.3-in. diameter IM-7 COPV was subjected to an ILH pressure 

schedule at LR « 0.3 to 0.9 




Pressure & Events vs. Time 

0 to 1 7500 s 1 7500 to 37500 s <“"• ) 



Event No. (filtered data) 


Results & Discussion 

COPVs 



IM-7 COPV data (^symbols) overlap IM-7 tow data 

(open green symbols) 



unloaded 


rupture 


Least squares fits (solid lines) and 99 % confidence intervals 
(dash-dot-dot lines) also shown for T 1 000 and Kevlar® 49 
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Results & Discussion 

COPVs 


• Lifetime (burst pressure) prediction of a COPY using a 
semi-empirical FR approach § 


Subject C/Ep connposite to 
internnittent load hold 
pressure schedule (red) 



Measure cumulative composite damage by 
monitoring approach of FR\o FR* 


First 8 Data Points: 



1.000 ± 0.010 

'v 

= 7,5401 75 ps 

m a^«ement ivuth hydroburst data (below) 

r ^ 

7,529 psi 
based on hydroburst data 
on2identKal COPVs 

All Data: 

Pfi’«0.961 +0.018 

f*V,.,.«.« = 7,870 ± 144 psi 
higher than imtially projected 7,S40ps value 


7,869 psi 
hi^erthan 7,S29psi value 
based on hydroburst data 
on 2 identKal COPVs 


1. Nichols, C., J. Waller, and R. Saulsberry, “Acoustic Emission Lifetime Estimation for Carbon Fiber/Epoxy Composite Overwrapped Pressure 

Vessels,” USRP Final Report, NASA-JSC Whites Sands Test Facility, Las Cruces, NM, August 2010. 

2. Waller, J., C. T. Nichols, D. J. Wentzel, and R. L. Saulsberry, Use of Modal Acoustic Emission to Monitor Damage Progression in Carbon 

Fiber/Epoxy Composites, QNDE Conference, July 18-23, 2010, San Diego, CA. 
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Results & Discussion 

COPVs 


AE event decay rate analysis on load holds using ASME Section X, Appendix 8 § 


y = Ae‘"+C 

Acceptance criteria from ASME 
Section X Mandatory Appendix 8 and 
NBlO-0601 Supplements : 

Acceptable Event Stability: 

- 0.1 <5 <- 0.0001 & R 2 > 0.80 


Observed acceptance criteria in WSTF 
IR&D IM7 COPV tests (more stringent): 

Acceptable Event Stability: 

-0.0030 < i? < -0.0019 & R2>0.90 


Shape factor B can also be expressed 
os the time required for the structure 
to emit 99% of events on a dwell. 

In(O.Ol) 25 to 40 minutes 
^99% ~ ^ (1535 to 2424 sec) 


§ 


Stability Analysis for ISS COPV S/N 008 

(Ch. 9) 


o Cum. Events at Dwell Observed UCL, B = -0.0019 

Observed LCL, B = -0.0030 Pressure 

Passing Test, B = -0.002209 



-600 0 600 1200 1800 


Seconds after the hold is reached (10 min. intervals) 


ASME Boiler and Pressure Vessel Code, Section X: Fiber-Reinforced Plastic Pressure Vessels, Section X, Appendix 8-620 Supplementary Examination 
Requirements (latest revision). 
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Conclusions 

COPVs 


• ASTM-based ILH methods were found to give a reproducible, 
quantitative estimate of the stress threshold at which significant 
accumulated damage began to occur. 

- FR AE events are low energy (<2 V^ias) 

- FR AE events occur close to the observed failure locus 

- FR AE events consist of more than 30% fiber breakage (>300 kHz) 

- FR AE events show a consistent hierarchy of cooperative damage for 
composite tow, and for the COPV tested, regardless of applied load 

• Application of ILH or related stress profiles could lead to robust 
pass/fail acceptance criteria based on the FR 

• Initial application of FR and FFT analysis of AE data acquired on 
COPVs is promising 

• EWMA knee methods for determining the ‘knee’ look very promising 

• Exponential AE event decay rate analysis on load holds look very 
promising 
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Back-Ups 


Load, Ibf 


Felicity ratio (FR) 



FR = 


stress at onset of significant acoustic emission during loading 

maximum previous stress plateau 
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Acoustic Emission Testing 


Acoustic Emission refers to the generation of transient elastic 
waves produced by a sudden redistribution of stress in a material. 
When a structure is subjected to an external stimulus (change in 
pressure, load, or temperature), localized sources trigger the 
release of energy, in the form of stress waves, which propagate 
to the surface and are recorded by sensors. 

(http ://www. ndt-ed . o rg/) 



•100 mV 
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Felicity Ratio 


Results & Discussion 



• For Kevlar-epoxy 4650 denier tow, correlation 

coefficients for ILH1 & 2 methods indicated good (R^ = 
0.866) to excellent (R^ = 0.985) agreement: 


200 250 

Max Prev. Load, Ibf 


Kevlar® 49-epoxy 



y = -0.0033x + 1.782 
R2 = 0.9048 


♦ Felicity Ratio 
Linear (Felicity Ratio) 



150 



200 


250 


y = -0.0023x + 1.545 
R^ = 0.985 


♦ Felicity Ratio 
— Linear (Felicity Ratio) 


300 


y = -0.000x+ 1.198 
R^ = 0.866 

♦ Felicity Ratio 
Linear (Felicity Ratio) 
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Results Sk Discussion 


• Characteristics of significant AE: 

- For Kevlar-epoxy, and T1000 and IM-7 carbon-epoxy, 
nonlinear increases in AE event rate were observed 
immediately before rupture, indicative of ‘critically intense’ AE 




Areas of critically intense AE activity also showed greatest 
violation of Kaiser effect, hence, the lowest FR values 
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Results & Di 


• IM-7 low vs. high energy events 


E < 2 V^-|^sec 


E > 2 V2-|asec 


• Low energy events have similar damage ‘footprint’ (top), while high 
energy events have a more variable damage ‘footprint’ (bottom) 

• Similar observation of a of a fiber breakage dominated “footprint” for 
FR events 
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Results & Discussion 

composite tow 


• Response surface generated using rotating mean method to identify where 
during in damage cycle best prediction of catastrophic failure can be made: 
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Results & Discussion 

COPVs 



WSTFArea 270 Test, Week 1, T1000 COPVs (General Dynamics) 
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